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We study from the point of view of quantum logic the 
properties of the collective oscillations of a linear chain of 
ions trapped in a linear Paul trap and composed of two ion 
species. We discuss extensively sympathetic cooling of the 
chain and the effect of anharmonicity on laser-cooling and 
quantum-information processing. 

PACS: 03.67.Lx,32.80.Pj,42.50.Vk 



I. INTRODUCTION 

The rapid development of trapping techniques for neu- 
tral and charged particles has constituted a breakthrough 
in the investigation of quantum mechanical systems |Q. 
Among the many interesting experiments, ordered struc- 
tures of charged ions have been achieved in Paul and Pen- 
ning traps P|J. Such structures are composed from few 
up to thousands of particles, and they originate at low 
temperature from the combined effect of the Coulomb 
repulsive interaction among the ions and the trapping 
potential Q]. Therefore, their geometry depends intrin- 
sically on the trap set-up. 

The field of interest of these so-called " Coulomb Crys- 
tals" is rather broad, and in quantum optics they find an 
application in the ion-trap quantum computer 0. Here, 
a string of ions is proposed as a system for processing in- 
formation, using two stable or metastable internal states 
of the ions for storing the quantum information and co- 
herent interaction of the internal degrees of freedom of 
the single ions with the laser light for generating the uni- 
tary operations which process the information, while the 
coupling among ions is provided by the collective vibra- 
tional excitations of the chain. Present schemes for quan- 
tum information processing are based on the harmonicity 
of the ionic motion H|?J . This regime can be achieved by 
laser cooling || the string of ions, and to a good ex- 
tent the ions can be considered to vibrate harmonically 
around their equilibrium positions. 

The coupling to the environment gives rise to phenomena 
which destroy the quantum coherence required for infor- 
mation processing: Decoherence affects both the internal 
and the motional quantum states. For motional states, 
decoherence can be inhibited by applying laser-cooling 
to the ion-chain on a regular rate or even continuously. 
Laser-cooling does not destroy the quantum information 
stored in the internal states, provided only some ions are 
addressed by the cooling laser {cooling ions), while the 



chain is sympathetically cooled via the Coulomb interac- 
tion. Furthermore it allows for simultaneous information 
processing on the other ions (qu-bits), provided the quan- 
tum gates do not require quantum coherence among the 
vibrational levels, as for the gate proposed in pj and re- 
alized in ||. In this respect, it is rather difficult to find 
a candidate ion which is, at the same time, a good qu- 
bit and cooling ion. In addition, some realizations of ion 
strings for quantum information do not rely on the spa- 
tial resolution of the ions with the laser g|. Hence, one 
of the most recent issues in the ion-trap quantum com- 
puter is to use two different ionic species which compose 
the ion chain, one for quantum information processing, 
the other for laser-cooling jlO],[ll]]. This type of crystal, 
which we call here the two-species Coulomb crystal, has 
been already used for imaging the mechanical effect of 
radiation pressure on the crystallized ions and for sym- 
pathetic cooling of big ordered ionic structures . 
In this work we study the mechanical motion in a two- 
species linear crystal from the point of view of quantum 
logic, illustrating the general features and the differences 
from a linear crystal composed of ions of equal masses, 
and we discuss in some detail sympathetic cooling of the 
chain. The issue of decoherence is also discussed in con- 
nection to the mechanical properties of the system. In 
particular, in this paper we consider the decoherence due 
to the coupling of the ionic motion to the fluctuations of 
the electric field (which we assume to have an instanta- 
neous value with zero spatial gradient) along the crystal, 
and thus couples with the center of mass motion. Finally, 
we discuss the harmonic approximation and the effect of 
the anharmonic corrections on laser cooling and quantum 
information processing. An analogous analysis has been 
presented in pd| ] , where they studied one specific scheme 
for sympathetic cooling. Here, we study the more gen- 
eral system and its properties, and seek schemes which 
are more suitable for realizing quantum logic. In partic- 
ular, we provide some examples calculated with 115 In + 
and 25 Mg+ ions, which have been currently trapped and 
cooled in Garching ]l^ , |T4| ]. 

The paper is organized as follows. In Section II the small 
oscillations formalism is introduced, and the mechanical 
properties of a chain of ions are discussed in detail. In 
particular, sympathetic cooling of the two-species linear 
chain is studied, and the rates of cooling for different crys- 
tal configurations are derived. In Section III the effect of 
the anharmonic corrections on cooling and quantum logic 
is discussed. Summary and conclusions are presented in 
Section IV. 
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II. SMALL OSCILLATIONS 

In this section we investigate the collective vibrations 
around the equilibrium position of a chain of ions con- 
fined in a linear Paul trap, following the lines of the lit- 
erature of small oscillations p5[ . We assume the motion 
to be onc-dimcnsional, i.e. confined along the trap axis 
of the linear Paul trap. This corresponds to assuming a 
very steep radial potential |l6| . 

We consider a one-dimensional string of N ions with 
charge e and mass equal either to M or rn, which are 
aligned along the x axis, which corresponds to the axis 
of the linear Paul trap. Indicating with i the position of 
the ion in the chain (i = 1, N), the sequence of ionic 
masses is described by the array m = (mi, ...,tojv) with 
mi = M, m. The ions are confined by the electrostatic 
potential V$ and interact via the Coulomb repulsion [jl7| . 
Sufficiently far away from the electrodes, Vs can be con- 
sidered harmonic and the total potential V has the form: 
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where Xi is the coordinate of the ion i and uq is a constant 
with the dimensions of an energy over a distance squared. 
If the ions are sufficiently cold M , they crystallize around 

the classical equilibrium positions xf ^ , which are the so- 
lutions of the set of equations dV/dxi\ x ( ) — 0. Those 

solutions are independent of the mass, as the potential 
of the electrodes interacts only with the ionic charges. 
A characteristic quantity is the equilibrium distance be- 
tween two ions: 
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where the ions are displaced symmetrically with respect 
to the center of the trap. This quantity scales the inter- 
ionic distance for a N ion chain @. 

Assuming that the collective motion of the ions around 
the equilibrium position is harmonic, we approximate V 
with its Taylor expansion around xf* truncated to sec- 
ond order. The dynamics of the system are described by 
the Lagrangian 
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where qi — Xi—x^ are the displacements of the ions from 
the equilibrium positions, and Vij are real coefficients 
which have the form: 
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From (|3|) the equations for the normal modes of the mo- 
tion are 
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Vijp? = Kmtf? with a = 1, N 



(5) 



where the eigenvalues X a are real given the hermiticity of 
V^, and where (3 a is eigenvector at A Q . Stable and har- 
monic oscillations exist if the condition A Q > is fulfilled 
for any a, as it occurs in this case. Under this condition 
the frequency of the normal mode fl a is Q a = \f\Z. The 
eigenvectors f3f are orthogonal in the Riemannian met- 
ric with metric tensor M, where M is a diagonal matrix 
whose diagonal corresponds to the array m. For X a > 0, 
introducing the mass-weighted coordinates q[ — ^Jmlqi 
the eigenvalue problem can be rewritten as 
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where now V{j — Vij / \jmimj, and the metric tensor is 
the identity matrix, as for cartesian orthogonal coordi- 
nates. The eigenvalue problem is now equivalent to the 
one of N identical ions of unitary mass. The matrix {j3f } 
defines an orthogonal transformation, which reduces the 
system to the principal axes 7r Q of Vij and of the ki- 
netic term: 7r Q = J2i ftt'li- ^ n this representation the 
Lagrangian describes a set of N independent harmonic 
oscillators with frequencies Q a . We quantize the motion 
by associating a quantum mechanical oscillator with each 
mode. Then, denoting a a , a^ a the annihilation, creation 
operators for the mode a, respectively, the coordinate Tr a 
associated with the oscillator of frequency fl a is written 
as 7r Q = y/h/2n a (a a +«4j- Going back to the origi- 
nal set of oblique coordinates , they have the quantized 
form 
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Some general features can now be recognized. From 
Eq. (jTJ) it is evident that the eigenmodes depend on the 
values of the ionic masses. Furthermore, since the matrix 
V.j is symmetrical by exchange of any pair of ions, the 
properties of the motion will be mainly characterized by 
the symmetries of the sequence m. These properties are 
reflected in the eigenmodes of the motion {/3 Q }, and thus 
they affect the coupling of the crystal to radiation. 
We discuss these points below. First, we consider the 
properties connected to two different values of the ionic 
masses, by analysing the case of two ions. Then, we 
discuss the ones connected to the symmetries of m by 
considering a three-ion crystal. Finally, on the basis of 
Eq. (0) we study the mechanical effect of radiation on 
the crystal, and in particular sympathetic cooling of the 
chain. 



2 



A. Two ions of different masses 

We analyse here the two-ion crystal, where the ions 
have masses m and M — fim with fi real parameter, fi > 
1. For this case the secular equation @ can be solved 
analytically and the eigcnfrcqucncics of the motion have 
the form: 



to \ fl V A* 
with corresponding displacements: 
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where the first and second components refer to the par- 
ticles of mass m and M, respectively. Here, N± are 
the normalization factor, according to the scalar prod- 
uct ^(1)^ (1) + fiq i (2)q j (2) = Sij with i, j = ±. 
In Fig. 1 (a) and (b) we plot the eigenfrequencies and the 
eigenvectors q±, respectively, as a function of fi. The ra- 
tio fi = 1 corresponds to the the well-known case of two 
ions of equal masses in a linear trap, where the ratio of 
the eigenfrequencies are in the relation 1 : v3- For this 
value of fi, f2_ and f2+ correspond to the center of mass 
(COM) and stretch mode frequencies, respectively, as it 
can also be verified from Eq. (^). As fi increases, the 
value of the eigenfrequencies decreases and tends asymp- 
totically to the values — > and Q + — > ^2u Q /m. The 
limit r2_ — > corresponds to the case where both ions 
stand still at their equilibrium position, as it can be seen 
in Fig. 1(b), while for the limit + — > \j2u$lm the heavy 
ion does not move, and the light ion oscillates around its 
equilibrium position. 

In the following we will concentrate on the case fi > 1 . As 
it is apparent from Fig. 1 (b) and Eq. (0) the two modes 
preserve some characteristics of the case of two ions with 
equal masses: In the mode of eigenfrequency J7_ the ions 
oscillate in phase, whereas for i7 + they oscillate with op- 
posite phases. The two modes, however, do not corre- 
spond to the COM and relative motion any longer. This 
can be understood by observing that, in absence of in- 
teractions, the trap frequency for an ion of mass m is 
proportional to 1/ \/m. This argument applies to an TV- 
ion chain of two (or more) species, and can be verified by 
substituting the vector g COM = (1, 1, l)/y/N describ- 
ing the center of mass motion inside the secular equation 
(§); one obtains £\ = u = mi X COM for i = 1, N, 
which cannot be fulfilled for any value of A , unless 
all masses rrii are equal. 

The non-separability of the modes into center of mass 
and relative motion has some consequences on the dy- 
namics. For example, for two ions the anharmonicity 
(i.e. the corrections to the harmonic approximation of 
the potential in Eq. (||)) couples the two modes, whereas 
in the crystals of ions with equal masses the COM mo- 
tion is an exact eigenmode of the problem. A further 



consequence is the coupling of both modes to the fluc- 
tuations of the electric field at the trap-electrodes, since 
none of the modes is orthogonal to the COM motion. The 
strength with which each mode couples to this source of 
decoherence is a function of the mass-ratio between the 
species fi, as has been discussed in H]. 



B. iV-ion crystal 

As discussed above, the characteristic properties of the 
motion of a two-ion crystal are a function of the mass ra- 
tio fi. For crystals with N > 2 ions, some further param- 
eters characterize the properties of the motion: The num- 
ber of ions of each species and the sequence in which they 
are arranged. These two features are described through 
the array m. In one dimension, the relevant symmetry 
property of the sequence m is the symmetry under reflec- 
tions with respect to the center of the trap (which is also 
the center of the string). This corresponds to an invari- 
ance of the Hamiltonian under spatial parity transforma- 
tions. Be n> N ' the parity operator, defined on the wave 
functions \4>(x\, X2, xm)) of the N-ion Hilbert space as 
U^\(j>(x 1 ,X2,...,x N )} = \<f>(-Xi,-X2,—,-x N )). This 
operator has eigenvalues p = +1 (even), p = — 1 (odd), 
corresponding to the states with even and odd parity, re- 
spectively. If the array m is symmetric under reflections, 
the Hamiltonian for the small oscillations commute with 
n( w ) and the eigenmodes of Eq. (JsJ) are also eigenvectors 
of U( N > at the eigenvalue either p — 1 or p = —1. From 
a simple evaluation of the number of degrees of freedom, 
one can verify that the even modes are N/2 for an even 
number N of ions and (N—l)/2 for N odd. In particular, 
for symmetry reasons the central ion does not move in 
the modes of even parity of a chain with an odd number 
of ions. Thus, these modes are independent of the mass 
of the central ion, as it can be deduced from Eq. (||). 

It is instructive to take a closer look at the eigenfre- 
quencies of a chain of N ions as a function of all possible 
sequences m. We discuss the case of 3 ions, since it shares 
some similarities with the normal modes of a triatomic 
molecule, as discussed in textbooks [ p!5[ , and it exhibits 
features which can be extended to chains of large N. 
In Fig. 2 we plot the eigenfrequencies of a N = 3 chain as 
a function of all possible sequences of Indium and Mag- 
nesium ions, where the sequences have been ordered with 
increasing total mass Mc of the crystal. The mode of fre- 
quency r2i (solid line) is characterized by the oscillation 
in phase of the three ions. For ions of equal mass it cor- 
responds to the center of mass mode, but in all cases still 
represents the mechanical response of the whole crystal 
to excitations: In fact, f2i oc I / \J Mc and in general it 
does not show an appreciable dependence on the order 
in which the ions are arranged. 

The properties of the higher excitations depend on fi and 
on the sequence. In particular, the distance among the 
eigenenergies changes depending on where the heavy ions 
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are placed in the sequence. In addition, symmetric se- 
quences preserve the properties of a chain of three ions 
of equal masses. Thus, the eigenmode of frequency 0.2 
(dashed line) is characterized by the out-of-phase oscilla- 
tion of the external ions, whereas the central ion stands 
still. Hence, £^2 takes the same value for the sequences 
A and C, and for the sequences D and F. In asymmet- 
ric sequences the external ions still oscillate out-of-phase, 
but the oscillation amplitude of the central ion is large, 
independently of its mass, provided that // ^ 1. 
In the symmetric sequences of the eigenmode Q3 the ex- 
ternal ions oscillate in phase, whereas the central ion is 
out-of-phase, and its amplitude is a monotonic function 
of l//i. For the sequence C, in particular, < ^3. The 
two frequencies are almost degenerate since both modes 
correspond to the case where the outer ions move sym- 
metrically with respect to the center, since the central 
ion in ^3 has small displacements. In asymmetric se- 
quences (B,E) the light ions have large oscillation ampli- 
tudes, whereas the displacements of the heavy ions are 
smaller. 

These properties have some immediate implications for 
quantum logic with a two-component chain. For ex- 
ample, the eigenmodes of even parity of symmetric se- 
quences are decoupled from the fluctuations of the elec- 
tric field at the electrodes, and thus are good candidate 
for the quantum bus. This issue have been discussed 
quantitatively for a particular sequence in fll"|| . 
Furthermore, in symmetric sequences the parity opera- 
tor commutes with the one-dimensional Hamilto- 
nian H: Hence, the eigenmodes with odd parity are not 
coupled via anharmonicity to the eigenmodes with even 
parity. The axial motion, however, is coupled to the ra- 
dial motion by non-linearities, and in three-dimension 
there are no groups of modes which are decoupled from 
the others. We discuss this point in section III. 
Another important implication regards the spacing be- 
tween the eigenfrequencies. In fact, when choosing a 
mode for quantum logic, the distance in energy among 
the eigenfrequencies should be taken into account, since 
the presence of quasi-degeneracies will lower the effi- 
ciency of single mode addressing in quantum logic op- 
erations. Thus, it is not convenient to use the mode with 
eigenfrequency in the sequence C as quantum bus, 
given its closeness to fo- 
rmally, from the mechanical properties of the chain we 
can infer in which position the cooling ions should be 
placed for optimal sympathetic cooling of the chain. We 
analyse this aspect in the next subsection. 



C. Interaction of the crystal with light 

In the dipole approximation the coupling of the ex- 
ternal degrees of freedom of an atom with radiation is 
represented by the kick operator exp(ik • r), where k is 
the wave vector of light and r the atomic position. In a 



Coulomb crystal of atomic ions, optical light couples to 
the internal degrees of freedom of a single ion, which we 
consider here a two-level dipole transition, and to the ex- 
ternal degrees of freedom of the collective motion. Thus, 
assuming that the ion j scatters a photon and that \ipi) 
is the motional state of the crystal before the scattering, 
after the scattering the motion of the crystal is described 
by the state \ipf) given by: 



\ipf) = exp(ikxj)\i/ji) 
= e 1 ^ exp(ikqj)\ipi 

M 



(10) 



where </> = kx)j is a real scalar, and k is the projection 
of k on the axis of the crystal. The coupling of radiation 
to the collective modes of the crystal is visible by substi- 
tuting Eq. (0) into ([H]). Thus, the kick operator can be 
written as 



exp(ifc^) = n Q=1 e 



(11) 



where 77" is the Lamb-Dicke parameter for the mode a 
and the ion j, and is defined as: 
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Through the Lamb-Dickc parameter we can infer the me- 
chanical response of the crystal to the scattering of a 
photon. For one ion of mass m, it corresponds to the 
square root of the ratio between the recoil frequency 
lor = hk 2 /2m and the trap frequency fl: rj = y/ujn/Vl, 
and it determines how many motional levels are coupled 
by the scattering of one photon. The Lamb-Dickc regime 
corresponds to the case in which ur <C f2, and math- 
ematically to the condition r\^fn <C 1, where n is the 
vibrational state number. In this limit the kick opera- 
tor can be expanded in powers of 77, and a change in the 
quantum motional state due to the incoherent scattering 
of one photon is of higher order in 77 19 1. In this regime 
an ion can be sideband cooled to the ground state of the 
vibrational motion f2(i|j . 

In presence of more than one ion, the Lamb-Dicke pa- 
rameter 77" describes how the displacement of the ion j 
couples to the mode a. In particular, it determines (i) 
the possibility of addressing a single motional sideband, 
which appears when scanning a probe beam through the 
resonance frequency (and thus of exciting one mode se- 
lectively) HJ], and (ii) the possibility of laser-cooling a 
mode to its vibrational ground state, in analogy to the 
one-ion case. 

Let us first consider the response to light of a two-ion 
crystal, and compare the case where the ions have equal 
masses, as discussed in pl| , with the case where they 
have different masses. In the first case, due to the sym- 
metry of the configuration the Lamb-Dicke parameters 
for each ion are equal (apart for some difference in the 
sign). In the second case, we can see in Eq. ( |l2| ) that 
the Lamb-Dickc parameters depend on the mass of the 
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ion, and given the asymmetry of the crystal the geomet- 
rical factor has different values for the two ions. One 
might be tempted to think that the Lamb-Dicke regime 
can be more easily accessed by addressing the heavier 
ions. This, however, is actually not true for all modes. 
This consideration is particularly applicable to the mode 
of lowest frequency, which shares some properties with 
the center of mass motion, and in general describes the 
response of the crystal as a whole to the mechanical ex- 
citation. For this mode, the displacements of the two 
ions are comparable, and actually the displacement of 
the heavier ion is slightly larger, as can be seen in Fig. 
1 (b) . On the contrary, the eigenmode with frequency f2+ 
is characterized by smaller displacements of the heavy 
ion than of the lighter one. For this mode, the heavy ion 
might be well in the Lamb-Dicke regime, while the light 
ion is not. 

This situation can be visualized by comparing the ab- 
sorption spectra obtained by shining light on each ion 
separately. We define the absorption spectrum in a two- 
level transition with resonant frequency ojq and driven 
by a laser of frequency ljl through the function 1(5), 
where 5 = ojl — loq is the detuning, which is evaluated by 
summing all contributions to laser-excited transitions at 
frequency lul Ejj]: 



1(5) 



|(n|exp(z^)|l)| 2 P(n) 



(13) 



the driven ion couples to the mode to cool. It thus con- 
tains some information on where the cooling ions should 
be placed in the sequence so as to achieve more efficient 
cooling. This can be theoretically illustrated by a rate 
equation, which describes cooling of one mode in a chain 
in the Lamb-Dicke regime [p2[ . Here, we define the Lamb- 
Dicke regime with the condition: 



a|2„ 



<1 for a = 1, 



(14) 



where n a denotes the occupation of the mode a and 
where the set of ions {j} represents the set of positions 
of the cooling ions in the chain. Assuming that the laser 
interacts with the internal two-level transition of the in- 
dividual atom, in second order perturbation theory in the 
parameter 5/7, with g Rabi frequency and 7 decay rate 
of the atomic transition, the excited state can be elim- 
inated, and one obtains a set of equations projected on 
the electronic ground state, where populations and co- 
herences between different motional states are coupled. 
In the Lamb-Dicke limit such coupling can be neglected, 
thus reducing the equation for the TV-ion density matrix 
p in the low saturation limit to rate equations pi| . Fur- 
thermore, introducing the reduced density matrix p a for 
the mode a defined by: 



Pa 
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Here |n) = | (711,77,2)) are the motional states of energy 
E n = nihQi + n-ihSlii with n a occupation of the mode 
Q a (a = 1, 2), and P(n) is a normalized distribution over 
the motional states |n). The coordinate qj characterizes 
the driven ion. In Fig. 3 we plot 1(5) as a function 
of 5 for a crystal composed of a Magnesium and an In- 
dium ion. Figure 3 shows the absorption spectrum when 
driving separately the Magnesium and the Indium ion, 
for a thermal distribution over the motional states with 
average total energy E — 5hQ-, where for the chosen 
parameters fl- = 0.552 MHz and Q + = 1.456 MHz. In 
both cases the motional sidebands of the mode with fre- 
quency f2_ are visible, whereas when driving the In + ion 
the sidebands of Q+ almost disappear. In fact, as eval- 
uated from ( |l2|) the Lamb-Dicke parameters for are 
7 ?M g = 0.22, r]^ n = 0.38, whereas the ones for Q + are 
^Mg = 0-5, 77^ = —0.06, i.e. the weight of the motional 
sidebands for the mode Q + in Indium is two orders of 
magnitude smaller than the ones for the mode 0_. It 
is interesting to compare this result with the case of two 
ions of equal masses. In that case the absorption spec- 
trum is the same independently of which ion of the chain 
is driven. Then, if the COM mode is in the Lamb-Dicke 
regime, the relative motion mode is also in the Lamb- 
Dicke regime, since in that case the Lamb-Dicke parame- 
ter scales simply as the inverse of the squared root of the 
eigenfrequency. 

For crystals with N > 2 ions the factor /3" in ( |l2| ) re- 
flects the structure of the chain and consequentely how 



with /3i, (3 n -i 7^ ct, one can derive the rate equation 
for cooling of the mode a in one-dimension p3| : 



^-P(n a ) = g -Y,*)f [(na + l)A°LP(n a + 1) 



dt ■ 



(16) 
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- ((n a + l)A% + n a A a _)P(n a ) + n a A%P(n a - 1) 

the coupling with the other modes being of higher or- 
der in the Lamb-Dicke parameter. Here, P(n a ) — 
(n a I p a I n a ) and the coefficients A± are defined as 



Al = 
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(17) 
(18) 



where we have assumed that the Rabi coupling is spa- 
tially constant over the whole crystal and equal to g, 
and that the laser is tuned on the first red sideband of 
the mode !7 Q . Eq. ([l6|) is the sum of all contributions to 
cooling of mode a from the coupling g of the cooling laser 
to the driven ions. It is fully equivalent to the equation 
for cooling of one ion in a harmonic trap of frequency fl a 
in the Lamb-Dicke regime jlj|, apart from the scaling 
factor multiplying the term on the RHS of (|l6|): 



E 



(19) 
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Provided that ( |l9| ) is different from zero, it does not af- 
fect the steady state, but scales the rate of cooling of 
the mode a. The factor W a represents the contribu- 
tion of the cooling ions in the array to the speed of the 
process. Obviously, the largest cooling rate is achieved 
when all ions are driven by the cooling laser. In that 
case, W a has the form W™ ax = J2f=i \Vj\ 2 = u R /Q. a , 
where ojr is the recoil frequency of the single ion. The 
rate of cooling of each mode scales according to the rela- 
tion W™ ax = W™ ax fL a > /fl a , and it scales with the mass 
to of the cooling ion as M / " lax oc -Jl/m since uir oc 1/m 
and Q a oc 1/ \fm. 

For quantum logic we are interested in employing only 
some ions of the chain for cooling. We look thus for 
the best sequence and mode for cooling, where for "best 
sequence" we intend a compromise between the highest 
number of ions for quantum logic (i.e. the lowest number 
of cooling ions) and the best cooling rate. 
In Fig. 4 we plot the factor W a vs all possible sequences 
of N — 3 ions, made up of Magnesium and/or Indium 
ions, where the cooling ions are in (a) Indium and (b) 
Magnesium ions. Note that between the curves in (a) and 
(b) there is a scaling factor corresponding, as expected, 
to the squared root of the ratio of the ionic masses. The 
curves in Fig. 4 can be easily interpreted by considering 
the properties of the modes, as discussed in the previous 
subsection. Thus, as the mode of frequency Oi is char- 
acterized by an oscillation in phase of all ions, and does 
not strongly depend on the sequence, the cooling rate 
increases as the number of cooling ions increases. On 
the other hand, the mode with frequency is mainly 
characterized by the oscillation of opposite phase of the 
ions placed externally. Thus, large rates of cooling are 
achieved when the cooling ions are placed in the exter- 
nal positions. In particular, the cooling rates of the se- 
quences A and C, and of the sequences D and F are 
equal. In fact, these configurations are symmetric un- 
der reflection and the mode with eigenfrequency Q2 has 
even parity. Thus, the central ion does not contribute to 
cooling nor to quantum information processing with the 
mode. Finally, in symmetric configurations, the mode 
of frequency f2 3 is characterized by large oscillations of 
the central ions, and a relatively large rate of cooling is 
obtained by simply placing the cooling ion in the center, 
as can be noticed for sequence C in Fig. 4(a) and se- 
quence D in Fig. 4(b). In asymmetric sequences (B,E) 
the rate of cooling in (a) is small, whereas in (b) is large, 
as expected from the considerations made in the previous 
subsection. 

For the case here considered, and on the basis of the me- 
chanical properties only, the best sequence of cooling is 
E provided that the cooling ion has lighter mass (Fig. 
4(b)). 

In general, we can conclude that by preparing certain 
sequences one can have efficient sympathetic cooling of 
some modes using a relatively small number of cooling 
ions. This characteristic does not depend strongly on the 



mass-ratio \i between the qubit and the cooling ions. In 
a crystal with large total mass Mc, however, the Lamb- 
Dicke regime condition can be accessed more easily. In 
this respect, it would be better to use heavier ions for 
quantum logic. 

Finally, sequences with an even number N of ions are 
to be preferred over sequences with odd N, so that all 
positions in the chain contribute either to quantum logic 
or to the cooling process. 

III. EFFECTS OF THE ANHARMONICITY ON 
COOLING AND QUANTUM LOGIC 

The harmonic approximation of the mechanical poten- 
tial in (P relies on the assumption that V possesses strict 
local minima, around which the motion is well-localized. 
In that case, the higher orders of its Taylor expansion 
are a small correction. For two ions those terms have the 
form (for X2 > X\): 

OO OO 2 I A 

sv = j:6v^ = j:(-ir e -^[ q2 - qi r w 

where 8V^- n ' is the n-th order correction. The effect of 
these terms, the so-called anharmonicity, consists in caus- 
ing shifts to the motional energies, and mixing the eigen- 
states of the normal modes. Such mixing is in general a 
small correction to the eigenstates of the ideal case, but 
it may become particularly enhanced because of quasi 
degeneracies among the motional energies of the states. 
In fact, the density of motional states of an N-ion chain 
in the interval of energy [E, E + 5E] is approximately 
D(E) cx E N ~ X . Thus, as the number N of ions increases 
and/or for larger values of the motional energies, the dy- 
namics of quasi-degenerate states are definetely affected 
by the anharmonicity. Here, we discuss the effects of the 
departures from the ideal harmonic system, first on side- 
band cooling and then on the efficiency of quantum logic 
gates. 

In sideband cooling the laser addresses the motional 
sideband of the mode to be cooled. Thus, cooling will 
be efficient as long as the shift in energy caused by the 
anharmonicity is much smaller than the frequency of one 
phonon of the mode addressed. On the other hand, the 
mixing between the eigenstates will constitute a thermal- 
ization effect among the modes, and it will not constitute 
an obstacle to cooling as long as all modes are at suffi- 
ciently low temperatures, so that 8V is a small correction 
to the whole system. 

To obtain some estimates, we evaluate the order of mag- 
nitude of the shift to the energy in first order pertur- 
bation theory, and ask for which range of values of the 
vibrational numbers sideband cooling may still work. In 
pl[ it has been shown that in the perturbative regime 
SVtvSVW |3 . Thus: 
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x \ x J 

where ao a = y/Ti/miQ, a , n a vibrational number of the 
mode a, and mj mass of the lighter ion. In deriv- 
ing we have assumed n a /£l a > np/Qp (a, (3 = 
1,2, ...,N). Taking two ions, one Indium and one Mag- 
nesium, f^i = 1MHz and the mass of Magnesium 25 Mg + , 
then \ (SV)\/h w 6 x 10 _3 n^ 2 r2i, which implies that first 
order perturbation theory holds for n\ <C 80 (n 2 <C 30). 
In this limit, a laser tuned on the first sideband to the red 
of the mode fii cools the system to the ground state. We 
have verified this conjecture numerically: We have taken 
a two- ion chain composed of an Indium and a Magnesium 
ion, and considered sideband cooling in the Lamb-Dickc 
regime of one of the modes, comparing the case in which 
the mechanical potential is fully harmonic with the case 
where the third and fourth orders in the anharmonic ex- 
pansion have been included. We have not noticed any sig- 
nificant difference between the two cases, and the system 
is cooled efficiently to the ground state. In particular, no 
visible effect could be interpreted as due to the anhar- 
monic coupling between quasi-degenerate states. How- 
ever, according to the above estimates, in the interval 
of states of the numerical calculation the spectrum of 
energy is not very "dense" . Then, in order to verify 
numerically the effect of anharmonocities in presence of 
quasi-degeneracies we take a system with exact degenera- 
cies, and more specifically with two modes of frequencies 
f^! = il and £! 2 = 20. Here, in the Lamb-Dicke regime 
a laser sideband-cools the mode of frequency fii. We 
compare the harmonic with the anharmonic case, where 
here we simply substitute the chosen values fii, 2 in the 
quantized form of the displacements of Eq. (Q) . In Fig. 5 
the average occupation number for the modes (a) f2i and 
(b) SI2 is plotted as a function of the time. The dashed 
and solid lines correspond to the harmonic and anhar- 
monic case, respectively. In the harmonic case the mode 
ill is cooled and O2 is like "frozen", since it is coupled 
to radiation at higher orders in the Lamb-Dicke param- 
eter |2lj . In the anharmonic case the rate of cooling of 
the mode addressed is slowed down, whereas the mode 
2 is simultaneously cooled: The system is cooled as a 
whole, but on a relatively slower time-scale. Thus, the 
two modes thermalize on a time-scale which is faster than 
the cooling one. 

With respect to quantum information processing, 
present quantum logic schemes with ions are based on 
the harmonic properties of the motion. Thus, they are af- 
fected both by the shift in energy and the mixing induced 
by the anharmonicity. Quantum gates which require the 
preparation of the system in the ground state of the mo- 
tion will be very weakly affected, since in that part of 
the spectrum the anharmonic perturbation is extremely 
small, and there are no quasi-degeneracies. On the other 
hand, the perturbation may affect the efficiency of hot 
gates, since they operate on higher-lying motional states. 
Thus, the speed of hot gates must be faster than the rate 



of anharmonic coupling between quasi-degenerate states. 
We can define a time-scale TA nn for the anharmonic per- 
turbation TAnh ~ h/\(SV)\- The typical duration of a 
quantum gate must be shorter than TAnh- In presence 
of degenerate states which are coupled by three-phonon 
transitions (i.e. for tt a ~ 2£1^), from ( ^lj ) TAnh ~ 10 fisec 
given fl a = 1 MHz and n a ~ 10. This estimate, which is 
rather worrying if compared with the duration of a quan- 
tum gate @,|25| , reflects the worst case, which might occur 
for certain sequences, large number of ions and large ex- 
citations. Note that, if the degenerate states are coupled 
by a four-phonon transition (i.e. for fl a ~ 30^), then 
TAnh ~ 600 /isec. In general, we expect this problem to 
arise when working with a large number of ions and for 
high excitations. It could be minimized by choosing sym- 
metric sequences: in that case the "effective density" of 
motional states which are coupled by the anharmonicities 
of the axial potential will decrease, since only states with 
the same parity will be coupled to each other, whereas the 
coupling with the radial degrees of freedom is of higher 
order. 

As a general rule, however, one should avoid degeneracies 
among the radial and the axial frequencies. 

IV. CONCLUSIONS 

We have studied the small oscillations behaviour of a 
two-component linear crystal, with particular emphasis 
on the applications to sympathetic cooling and quantum 
logic with the ions, and have discussed the effect of an- 
harmonicity on the operations of the ion-trap quantum 
computer. We have seen that higher efficiency in quan- 
tum logic and sympathetic cooling are achieved by select- 
ing the right ionic sequences. That raises the issue of how 
to prepare the desired sequence of ions. A rigorous in- 
vestigation in this direction should take into account the 
full non-linear potential in the three-dimensional space 
and it is subject of on-going investigations. 
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ported in parts by the European Commission within 
the TMR-networks ERB-FMRX-CT96-0087 and ERB- 
FMRX-CT96-0077. 
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FIG. 1. (a) Eigenfrequencies of the two-ion crystal as a 
function of t he mas s ratio fi. The frequencies are rescaled by 
the factor ^Juo/m. (b) Corresponding displacements in the 
oblique axes plotted as a function of /i. Solid line: particle of 
mass M. Dashed line: particle of mass m. 
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FIG. 2. Eigenfrequencies £l a vs all possible configurations 
of a N = 3 ionic sequence made up of Indium (black circles) 
and/or Magnesium (white circles) ions. The eigenfrequencies 
are rescaled by the value of fii for a chain of Magnesium ions. 
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FIG. 3. Absorption spectrum 1(5) vs the detuning 8 for 
a crystal of an Indium and Magnesium ion obtained shining 
light on Magnesium and Indium, for a thermal distribution 
with average energy 5hfl-. Here, J?^ g = 0.22, r)^ = 0.5, 
T]j n = 0.38, r/+ = -0.06, as follows from @. 
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FIG. 4. Cooling rate W a vs all the possible configurations 
of a iV = 3 ionic sequence made up of Indium (black circles) 
and/or Magnesium (white circles) ions, where W a has been 
rescaled by the value of for a chain of Magnesium ions. 

Solid line: Eigenmode Qi; Dashed line: Eigenmode ^2; Grey 
line: Eigenmode Q3 [cif. Fig. 3]. 
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FIG. 5. Average occupation of the mode (a) fii = Q 
and (2) ^2 = 2f2 as a function of the time in unit of the 
optical pumping time r, for the harmonic (dashed line) and 
anharmonic case (solid line). Here, = 1 MHz, r/n ~ 0.1, 
?72n = 0.1/V5, 7 = O.lfi, 5 = U) L - w = -fi, 9 = O.Olfi, 
r = 7/<? 2 . Insets: Population of the modes (a) Q,i and (b) Q.2 
as a function of their vibrational number at time t = 150r for 
the anharmonic case. 
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